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Abstract. I analyze unpublished or partially published archival ROSAT data of the p Ophiuchi dark cloud. This 
set of seven overlapping ROSAT HRI pointings, composed of eight ~one-hour exposures, detects mainly the X-ray 
brightest T Tauri stars of this star-forming region. Only two HRI sources are new X-ray sources, and their optical 
counterparts are proposed as new Weak T Tauri star candidates. Meanwhile the ROSAT HRI caught during just 
one exposure a weak X-ray source (£=10; SNR— 4.1 a for Gaussian statistics) among a group of three embedded 
young stellar objects including two Class I protostars. Previous ROSAT PSPC, ASCA GIS observations, and as 
I argue here one Einstein IPC observation, have already detected an X-ray source in this area, but this higher 
angular resolution data show clearly that X-rays are emitted by the Class I protostar YLW16A. This is the 
second Class I protostar detected by the ROSAT HRI in this dark cloud. The determination of the intrinsic X- 
ray luminosity of this event, Lx[0.1-2.4 keV] = (9.4-450) x 10^" ergs~^, critically depends on the source absorption 
estimate. Improvements will be obtained only by the direct determination of this parameter from fitting of Chandra 
and XMM-Newton spectra. 
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1. Introduction 



tions ([Grosso et al. 1997 , hereafter paper I; Neuhausci 



As in Llie 80's Lliu EinsLnin X-iay ubsyivaLuiy disuuv- 
ered the high X-ray variability of low-mass pre-main se- 
quence stars, called T Tauri stars, the ROSAT satellite 
reported in the 90 's the first clues of the X-ray activity of 
youngqi sLellai ubJucLs ('^ 10''' ^i), called Class I piuLusLais 



fc Preibisch 1997). With the goal of preparing the up- 



coming study of these two star-forming regions with the 
new generation of X-ray satellites, Chandra and XMM- 
Newton, their ROSAT observations must be definitively 
exploited. This work was already done by Neuhauser & 



( Feigel &ou &, MouLmcilc 1999 ). Class I pioLustdiS ( Lada 



Preibisch (1997) on R CrA. I study here the unpublished 



1991) are composite objects including a central forming 
star surrounded by an accretion disk ~10-100AU in ra- 
dius, and embedded in an extended infalling remnant en- 
velope of gas and dust up to 



or partially published observations of the p Ophiuchi dark 
cloud to search for X-rays from Class I protostars. 



10"^ AU in size ( 3hu ct al. 2. ROSAT observations 



1987). Using the Position Sensitive Proportional Counter 



(PSPC), [Casanova et al. (1995) reported in the p Ophiuchi 
star- forming region (^145 pc; [de Zeeuw et al. 1999| ) the de- 
tection of seven X-ray sources coinciding with Class I pro- 
tostars, but possible faint non-protostellar counterparts 
were also found in the same error boxes. The ASCA satel- 
lite detected more clearly a group of five Class I proto- 
stars in the R CrA star-forming region (~130pc), but the 
positional uncertainties was quite large (~ 20"), and the 



I found in the ROSAT archive that ten overlapping point- 
ings were performed on the p Ophiuchi dark cloud: one 
with the PSPC and nine with the HRI (see Table | for 
the log of these observations). Among these HRI point- 
ings a sequence of six shor t exposures (~ hour) s preaded 



individual source only partially resolved (Koyama et al 



over one week, was used by [Damiani ct al. (1996)| to study 
the X-ray variability of the T Tauri stars SR9 and SR12A- 
B. One more HRI pointing remains without published re- 
sults. I thus analyzed these data (8 segments in all; see 



199£). In both cases ROSAT, with its High Resolution 



Imager (HRI) providing a much better spatial resolu- 
tion (FWHMwS"), has helped to confirm these detec- 
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Table |1|) with EXSAS ( [Zimmermann et al. 1997] ). 

Source detection was performed on each segment sep- 
arately with the standard command DETECT/SOURCES, 
which generates a local source detection by a sliding- 
window technique followed by a maximum likelihood test, 
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Table 1. List of the ROSAT observations of the p Oph dark cloud. Only the segments #1-8 are analyzed here. 



ROSAT 


Proposal name 


PI name 


Exp. 


Start 


End 


rv Tonnn 

"-••J ZUUU 


u Tonnn 

'-'J ZUUU 


Ref. 


Segment 


cLrchive 






fksl 


[vvmmddl 


[vVTTlTTlddl 


16'' 


-24° 




41 

TT 


200045p-0 


p Oph Core 


Montmerle 


12.9 


910305 


913010 


26'" 31:0 


31'48" 


(1) 




200045p-l 


p Oph Core 


Montmerle 


19.9 


910908 


910908 


26'" 3 1=0 


31'48" 


(1) 




201709h 


SR9, SR12 


Damiani 


3.3 


940829 


940829 


27'" 28=1 


31'48" 


(2) 


1 


201710h 


SR9, SR12 


Damiani 


5.2 


940831 


940831 


27'"28n 


31'48" 


(2) 


2 




oxly, orxlz 


Damiani 


7 


0/1 nam 




97111903 1 








201712h 


SR9, SR12 


Damiani 


5.1 


940903 


940903 


27'"28=1 


31'48" 


(2) 


4 


201713h 


SR9, SR12 


Damiani 


8.1 


940904 


940904 


27'" 28=1 


31'48" 


(2) 


5 


201714h 


SR9, SR12 


Damiani 


3.7 


940905 


940905 


27'"28=1 


31'48" 


(2) 


6 


201618h-l 


ROX20 


Zinnecker 


1.9 


940917 


940917 


27™ 14=0 


51'36" 




7 


201834h 


Oph Core F 


Montmerle 


12.6 


950309 


950314 


27^26=0 


40'48" 


(3, 4) 




201618h-2 


ROX20 


Zinnecker 


3.6 


950817 


950817 


27^14=0 


51'36" 




8 


201834h-l 


p Oph Core F 


Montmerle 


27.8 


950818 


950820 


27'"26=0 


40'48" 


(3, 4) 




201835h 


p Oph Core A 


Montmerle 


51.7 


950829 


950912 


26'"02!0 


23'24" 


(4) 




201834h-2 


p Oph Core F 


Montmerle 


37.5 


960907 


960911 


27'"26=0 


40'48" 


(4) 





Note: In the R OSAT archive "p" (resp . "h ") is for the PSPC (re sp. HRI) instrument. 

References: (1) [Casanova et al. (1995) ; (2) Damiani et al. (1996); (3) paper I; (4) Grosso et al. (2000) 



which compares the observed count distribution to a 
model of the point spread function and the local back- 
ground to discriminate sources from statistical Poissonian 
background fluctuations. The likelihood of existence, de- 
fined as £=— InPo (with Pq the probability of the null 
hypothesis that the observed distribution of counts is only 
due to a statistical background fluctuation), provides a 
maximum likelihood measure for the source detection. 
I accepted only detections with £>10 {SNR>A.1<t for 
Gaussian statistics) to reduce the number of spurious de- 
tections per field to ^0.7-~1.0. Identification of these X- 
ray sources was made by cross-correlation with published 
list of confirmed or suspected cloud members (e.g. Andre 
& Montmerle 1994), IR surveys (e.g. Barsony et al. 1997), 
and optical catalogue (Monet et al. 1996). To correct X-ray 



positions for boresight errors, I selected the X-ray sources 
not associated with protostars having a positional error 
lower than 2" both in a and 5, and compared their posi- 
tions with their IR counterparts in the 2MASS catalogue 



(second incremental release; Cutri et al. 2000); the mean 



offsets in a and S was then substracted to the X-ray posi- 
tions of all X-ray sources, and the residual dispersion was 
quadratically added to the positional errors. Table || gives 
the X-ray source list. 

These short exposures detected mainly the X-ray 
brightest T Tauri stars of this star-forming region. Only 
two HRI sources are new X-ray sources. They are located 



~25' south (-Ipc) of the core E/F ( |Loren et al. 1990D 
and identified with optical stars]]: they are reliable Weak 
T Tauri star candidates. Meanwhile a weak HRI source 
(£=10), detected just in segment 8, is associated with the 
IRAS source YLW16 - a group of two Class I protostars 
(YLW16A; and YLW16B, also called IRS46) and one em- 
bedded Classical T Tauri star (GY262). In order to dis- 
cuss the possible identification of this weak HRI source 

^ 11589441 (B=14.8, i?=12.6); 11614157 (B=17.9, i?=14.1). 



20 



30 



LYLW16B/IRS46 (lA (ROSAT VsPC) , 



-24 39 40 



; X R0XA-9B 
1 (ASCA GIS) 

16''27"30= 




YLW16C/GY262 (II)- 



LFAM35 (VLA)\ 

, "\YLW16A/IRS44 (I) 
Raw ROSAT HRI/ x \ ^. , „„o.n, ,,L, 
error box\ X M^'na' ROSAT HRI 



error box 



27= 



Fig. 1. X-ray detection of the Class I protostar YLW16A 
with the ROSAT HRI. The positions of the IR sources 
in this area are from the 2MASS catalogue (Cutri et al. 
2000; o-sMASS ^ 0.3"). The position of YLW16D, given 
in the Barsony et al. (1997) coordinate frame by Lucas 
& Roche (1998), was converted to the 2MASS coordinate 
frame applying the position offset found for YLW16A in 
2MASS compared to Barsony et al. (1997). The position 
of LFAM35, the VLA detection of YLW16A (Leous et al. 
1991), is also shown. The circle radii indicate with the 
one cr positional error of the different observations. The 
ROSAT HRI error circle after boresight error correction 
includes only YLW16A. 



with protostars. Fig. |l| visualizes the different members of 
YLW16 in the 2MASS reference frame, to be consistent 
with my previous astrometric correction. The 2MASS po- 
sition of YLW16A is within the 4"-radius error circle of 
the HRI source. Using a Monte-Carlo simulation, I found 
that the probability that this identification is only due to 
chance is 1.6%. This reliable identification strenghtens the 
significance level of this weak HRI detection. 
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Table 2. X-ray sources detected by the ROSAT HRI 
with a likehhood of existence >10 in the segments #1-8. 
Col. (1) gives the IR/optical counterpart name with : GY 
= Greene & Young (1989); IRS = Wilking et al. (1989); 
ROXs = Bouvier & Appenzeller (1992); SR = Struve & 
Rudkj0bing (1949); YLW = Young et al. (1986); numbers 
are from the PMM USNO-Al.O catalogue (Monet et al. 
1996; the designation prefix 0600-1 was here cut away). 
Col. (2) indicates the IR classification. The better X-ray 
positions of this segment set in Col. (3-5) are corrected 
from boresight errors. The segment number as defined in 
Table ^ is given in Col. (6). C, in Col. (7), is the likelihood 
of existence of the X-ray sources for each observation. The 
count rates are given in the ROSAT 0.1-2.4keV energy 
band in Col. (8). 
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III? 


26' 


"04° 5 


-24°57'50" 
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20.2 ±3.0 


SR24N 


II 


26' 


"58=3 


-24°45'32" 


3 


7 


10.8 


4.6 ± 1.7 














8 


13.3 


2.4 ±0.9 


GY194 


III 


27' 


"04=8 


-24°42'18" 


5 


7 


10.1 


3.9 ± 1.6 


ROXs20B 


III 


27' 


"15=1 


-24°51'40" 


3 


7 


18.5 


4.0 ± 1.5 














8 


38.0 


4.5 ± 1.2 


SR12A-B 


III 


27' 


"19! 7 


-24°41'39" 


2 
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143.2 


16.7 ±2.4 














2 


123.0 


11.3 ± 1.6 














3 


138.6 


24.2 ±3.2 
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128.2 


12.3 ± 1.7 
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2150.2 


85.4 ±3.4 
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13.0 ±2.0 
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11.7±2.7 
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19.0 ±2.4 


YLW16A 


I 


27' 


"28=0 


-24°39'38" 
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2.2 ±0.9 


1614157 


III? 


27' 


"32^5 


-25° 06' 16" 
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21.3 


12.1 ±3.1 
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20.8 


8.2 ± 1.9 
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27' 


"33=1 


-24°41'20" 
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11.8 


2.3 ±0.8 
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10.7 


2.0 ±0.7 
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27' 


"38=2 


-24°36'58" 
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1.8 ±0.8 


SR9 


II 
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"40^4 
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29.5 ±3.2 
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6 
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19.5 


4.8 ± 1.4 














5 


29.8 


4.6 ±0.9 














6 


42.1 


7.4 ± 1.6 
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II 
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6 
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As shown in Fig. ^ other instruments detected also 
an X-ray source in this area. The ROSAT PSPC source 
ROXRl-45 (^ asanova et al. 1995D was identified with 
YLW16. The ASCA CIS source R0XA-9B ( |Kamata et 
al. 19^) was only identified with YLW16B, whereas 
YLW16A is also in the ASCA CIS p ositional error box 
as mentioned by [Carkner et al. (1998); moreov er the new 
IR source discovered by Lucas fc Roche (1998) , YLW16D, 
is in the same error box. An ASCA follow-up observa- 
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Fig. 2. Einstein IPC detection of the YLW16 IR source 
group. This intensity contour map is an enlargement of a 
convolution of the counts distribution (pixel size binned 
to 30") of an archival Einstein IPC observation of the p 
Ophiuchi dark cloud (8 sep. 1979; 2.1 ks exposure) with a 
Gaussian profile (FWHM=1'). The image background is 
i3=2.5xlO~'*ctss~^ arcmin"^. The intensity scale is lin- 
ear, and was choosen so that the first and fourth contours 
match the Montmerle et al. (1983)'s contours n—\ and 
4, and thus represent levels 2 x 1.5^+"i3 above standard 
background. The Einstein IPC FWHM in this area, lo- 
cated ^ 20' away from the axis, is larger than 1'. The 
arrow heads point the 2MASS positions of the IR coun- 
terpart of Einstein IPC sources (the ROX number is given 
after the sign "="). Asterisks (resp. cross) show the Class I 
protostar (resp. T Tauri star) positions (see also Fig. |l]). 
An Einstein X-ray source is associated with the YLW16 
IR source group. 



tion, made 3.5 years after, did not detect again ROXA- 
9B (Tsuboi et al. 1997): this X-ray source is variable. I 
noted also in only one of the Einstein IPC observations of 
Montmerle et al. (1983; see observation 13749.2 Fig. 1.2) 
an X-ray emission around YLW16, blended with the 
Einstein X-ray sources associated with SRI 2 (R0X21) and 
GY292 (ROX25), which was not in the Montmerle et al.'s 
detection list. To solve this point, I took from HEASARC 
the archived screened photon event list^, and revisited it 
using the package XANADU/XIMAGE (also available at 
HEASARC). I selected only the 8 Sep. 1979 events and 
energy bins 3-13 (corresponding to 0.3-8.2 keV). From 
this event list, I studied the X-ray source standing near 
YLW16 with the interactive command sosta, which al- 
lows to estimate the background from a nearby area free 
of sources, and to tune the source box size to exclude 
events from neighbouring sources. I found for this source 
a signal to noise ratio ~3.5ct, and a corrected intensity 
S'~0.02ctss~^. Fig. |2| shows a contour map constructed 
from these data. I conclude that the Einstein IPC also 
detected X-rays from YLW16 (which was not yet discov- 
ered from IR at this time), but the Einstein IPC FWHM 



Sum of 13749.1 (8 March 1979), and 13749.2 (8 Sep. 1979). 
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(~1') is too large to find unambiguously the IR counter- 
part. By contrast to these previous detections with other 
X-ray instruments, we can see clearly in Fig. ^, thanks to 
the better angular resolution of the HRI, that the X-ray 
emission detected by the ROSAT HRI comes only from 
the Class I protostar YLW16A. 

3. X-ray luminosity of YLW16A 

YLW16A is the second Class I protostar detected by the 
HRI in the p Ophiuchi dark cloud after YLW15 (paper I). 
It is remarkable to detect it with only an one-hour HRI ex- 
posure, whereas no detection was obtained with a 20 times 
longer HRI exposure (see Table |l|; [Crosso et al. 200(]| ). The 
fact that YLW16A was not detected by the HRI only one 
day before, with an exposure 2 times longer, implies a vari- 
ation of its X-ray luminosity by at least a factor 2 in less 
than one day. This suggests that YLW16A was detected 
during a state of higher X-ray luminosity, probably due to 
an X-ray flare, usual in young stellar objects (see Feigelson 
& Montmerle 1999). However the source X-ray variability 
cannoti be 

with only ~ » cts detected during this short exposure. It is 
impossible to deduce from this HRI detection whether the 
previous X-ray observations by other instruments actually 
detected YLW16A, as this X-ray source is not constant. 

The determination of the intrinsic X-ray luminosities 
of such embedded young stellar object critically depends 
on the absorption of the X-ray photons by the gas along 
the line of sight, A'h, combination of both the circumstel- 
lar material and the interstellar medium. As the HRI has 
no spectral resolution, no direct information on A'h can be 
provided by the X-ray event. I thus derive iVfj from the 
source visual extinction due to the dust, Ay, ass uming 



provide an accurate value of A^h to constrain the X-ray 
luminosity observed by the HRI. 
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